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ABSTRACT 

Context. The Red MSX Source (RMS) survey is a multi-wavelength campaign of follow-up observations of a colour- selected sample 
of candidate massive young stellar objects (MYSOs) in the galactic plane. This survey is returning the largest well-selected sample of 
MYSOs to date, while identifying other dust contaminant sources with similar mid-infrared colours including a large number of new 
ultra-compact (UC) Hii regions. 

Aims. To measure the far-infrared (IR) flux, which lies near the peak of the spectral energy distribution (SED) of MYSOs and UCHii 
regions, so that, together with distance information, the luminosity of these sources can be obtained. 

Methods. Less than 50 % of RMS sources are associated with IRAS point sources with detections at 60 fim and 100 yum, though the 
vast majority are visible in Spitzer MIPSGAL or IRAS Galaxy Atlas (IGA) images. However, standard aperture photometry is not 
appropriate for these data due to crowding of sources and strong spatially variable far-IR background emission in the galactic plane. 
A new technique using a 2-dimensional fit to the background in an annulus around each source is therefore used to obtain far-IR 
photometry for young RMS sources. 

Results. Far-IR fluxes are obtained for a total of 1 1 13 RMS candidates identified as young sources. Of these 734 have flux measure- 
ments using IGA 60 /jm and 100 yum images and 724 using MIPSGAL 70 yum images, with 345 having measurements in both data 
sets. 

Key words. Stars: Formation - Stars: Pre-Main Sequence - Hii Regions - Infrared: General - Techniques: Photometric - Surveys 
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^ 1. Introduction IWvnn-Williamsl (1 19821) and iHenning etaTI (11984 represented 

the largest collection of MYSOs but, with only of order 50 

O Massive (M > 8 M^) stars dominate the regions they form in due ^^^^^^^ between them, this is a factor of 10 smaller th an ex- 

L: to their strong ionising radiation and their powerful outflows and p^^^^^ (iLumsden et all ^OOl . This led i Lumsden et all (12001) 

^ winds. However, understanding the early stages of their forma- -^^^^^^ ^ ^^^^^^ ^ 2000 candidate MYSOs from the 

^ tion and evolution remains difficult due to their rarity, the em- ^idcourse Space Experiment (MSX) point source catalogue 

H beddedness of the formation process and the short timescale over (gganetaL .1999. .2003.) using colour selection criteria. 
5t which it takes place. Of particular interest is the massive young 
stellar object (MYSO) phase, which is mid-infrared (IR) bright 

but radio quiet, where core hydrogen burning has probably be- The Red MSX Source (RMS) survey is a campaign of 

gun but an Hii region has yet to form. As such, MYSOs repre- follow-up observations designed to identify other dusty ob- 

sent one of the earliest phases of the life of a massive star where Jects which contaminate this sample of MYSO candidates (ul- 

it is truly massive while major accretion may still be ongoing, tra compact (UC) Hii regions, low-mass YSOs, Evolved stars, 

as evidenced by their p owerful bipolar molecular outflows (e.g. proto-planetary nebula ie (PNe) arid PNe) and gain informa- 

iBeuther et al.Ll2002allq) . tion about the sources (iHoare et am2005l : lMottram et al.L 120061; 

A major stumbling block to studying MYSOs, and by ex- lUrquhart et al.1 |2008b|). These include; ~ 1" resolution con- 
tension massive star formation in general, is the lack of large tinuuiii radio observations to identify regions of ionised emis- 
well-selected samples of sources from which to differentiate sion ("UrquharteLaU |2Q07al |2009^; ground-based ~ 1" reso- 
the global properties of all such sources from individual pe- lotion mid-IR imaging where Spitzer IRAC (3.6 //m, 4.5 //m, 
culiarities and variations. Until recendy, the relatively small 5-8 ^m and S.O ^m) imaging from tfaeGLIMPSE survey 
catalogues of primarily serendipitously discovered MYSOs by (IBeniamin et alj |2003|; [Churchwell et alj |2009D is unavailable 

to identify multiplicity in the ~ 18" MSX beam, obtain bet- 

* The full version of table [2 is only available in electronic form at ter astrometry and ensure that the survey does not bias against 

the CDS via anonymous ftptocdsarc.u-strasbg.fr (130.79.125.5) or via YSOs near H ii regions (e.g. Mottram et al., 2007); '^^CO molec- 

http://cdsweb.u-strasbg.fr/cgi-b in/qcat7J/A-i- A/| ular line observations in orde r to obtain kinematic distances 

** E-mail:joe@astro.ex.ac.uk (lUrquhart et al.L l2007bl l2008al) ; and near-IR spectroscopy (e.g. 
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Fig. 1. A comparison of an IRAS ISSA 60 //m image (top), 
an IGA 60 jjm image after 20 HIRES iterations (centre) and a 
Spitzer MIPSGAL 70 fim (bottom) image of the same region 
centi-ed on G015. 00+00.00. The relevant features of IGA and 
MIPSGAL data are discussed in ^22]and ^231 



IClarke et al.L l2006l) to distinguish between YSOs and evolved 
stars. 

In order to obtain the luminosities of the young massive stars 
identified by the RMS survey, photometry is required over as 
much of the spectral energy distribution (SED) as possible. In 
particular, the SEDs of YSOs and UCH ii regions peak in the far- 



IR between ~ 60 fim and ~ 120 jjm. It is therefore important that 
far-IR photometry is available for those candidates identified as 
young, i.e. associated with star formation activity. Despite cov- 
ering the majority of th e sky, the IRAS point source catalogue 
(lBeichmanetaUll988l) does not contain entries for ~ 28 % of 
young RMS sources, and a similar number have only upper limit 
detections at 60 yum and/or 100 jjm. In addition, the background 
emission at far-IR wavelengths is strong and spatially variable 
within the galactic plane for IRAS data which adds complexity 
to background subtraction. Simple aperture photometry, which 
estimates the background emission using the average within an 
annulus surrounding the source is therefore unlikely to provide 
accurate measurements. However, fitting the point spread func- 
tion (PSF) of the detector is also difficult due to the crowded 
nature of sources in the galactic plane and interaction between 
the extended background and the PSF. 

In this paper, we begin by discussing in ^the RMS sam- 
ple of massive young stars and the data used to obtain far-IR 
fluxes for them: HIRES reprocessed IRAS Galaxy Atlas (IGA 
[Cao et al.', '1997) 60 fim and 100 yum images and MIPSGAL 
(Carev et al., 2009) 70 jum images. In |3]we outline a new tech- 
nique for performing aperture photometry using a 2-D fit to the 
background within an annulus, rather than a simple average. The 
results of photometric measurements using this technique are 
presented in ^ and discussed in ^ We then summarise and 
reach our conclusions in ^ 



2. Data 

2. 1 . The RMS sample of massive young stars 

Using the RMS follow-up observations and data from the liter- 
ature (e.g. 2MASS, GLIMPSE), aU of which are available on 
the RMS database (http://www.ast.leeds.ac.uk/RMS/), each can- 
didate MYSO has been assigned one of several identifications 
(see Mottram, 2008, for more details). Young sources can be as- 
signed 'YSO', 'Hiiregion', 'Hii/YSO' and 'Young/Old Star'. 

Sources assigned the 'YSO' identification are radio-quiet 
mid-IR point sources associated with '^CO emission which do 
not look isolated in near-IR images. Where near-IR spectroscopy 
is available, thes e sources usually show H i emission, though not 
at Case-B ratios dBaker & Menzellll938l) . and can also show H2, 
He I, Fe 11 and [Fe 11] emission, as well as CO bandhead emission 
and/or absorption. 'Hii region's are either extended in ~1" res- 
olution mid-IR imaging and/or radio loud, with Case-B H i line 
ratios in the near-IR. 'H 11 / YSO' sources either have indications 
that both a YSO and an H 11 region are present within the MSX 
beam (18"), or have conflicting evidence as to their identification 
with insufficient information currently available to determine be- 
tween the 'YSO' and 'Hii region' designations. Finally sources 
identified as 'Young / Old Star's are unlikely to be H 11 regions 
or PNe as they are radio quiet and unresolved in the mid-IR, but 
have insufficient evidence to determine whether the source is a 
YSO or an evolved star As more information becomes available 
identifications are updated where needed, particularly from our 
near-IR spectroscopy and the luminosities derived using the data 
presented in this paper, which will be the subject of a forthcom- 
ing paper (Mottram et al., in prep.). 



2.2. The IRAS Galaxy Atlas 

The IRAS sky survey scanned most regions of the sky more than 
once, often at different scan angles, t hus additional sp atial infor- 
mation is available within the data. ICao et all (1 19971) therefore 
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undertook HIRES repr ocessing of the or i ginal IRAS Sky Survey 
Atlas (ISSA) images (IWheelock et all 11994 for the galactic 
plane in order to produce the IGA. These data have improved 
resolution (1.0' x 1.7' at 60 yum, 1.7' x 2.2' at 100 fim) with re- 
spect to both the original ISSA, as shown in figure [1] and the 
more re cent Improved Reprocessing of the IR AS Survey (IRIS) 
images (Mivi lle-Deschenes & Lagachel l2005h . which have res- 
olutions between 3 .5' and 5'. Though the IR IS data are bet- 
ter calibrated due to lMiville-Deschenes & Laga che (2005) using 
COBE/DIRBE data, the improved resolution of the IGA is cru- 
cial when considering regions in the galactic plane. 

Due to the shape of the IRAS detector pixels and the detector 
layout used, the IRAS sky survey had different resolution par- 
allel and perpendicular to the scanning direction of the satellite. 
This caused the point source function (PSF) of raw IRAS images 
to be elliptical, rather than circular - a feature which is retained 
in the IGA (see figure[Tll, though IRIS images are smoothed to a 
circular PSF with resolution closer to the larger dirn ension (i.e. 
cross-scan, see lMiville-Deschenes & Lagacheil2005h . The angle 
between the major axis of the PSF and the x-axis of the image 
also varies across the IGA images due to the changing scan di- 
rection of the IRAS satellite across the sky as it orbited the earth 
and the earth orbited the sun. 

As a part of HIRES processing, beam sample map images 
were produced which have the PSF placed at the same posi- 
tions on each tile on a smoothed version of the 20 iteration im- 
age in order to provide a reasonable estimate of the local back- 
ground. These were used to measure the local beam angle for 
each source. In addition, HIRES processing results in a neg- 
ative ringing artifact around the central maximum of the PSF 
(see figure [T]) which must be avoided during aperture photom- 
etry, caused by i nteraction betwe en the PSF and the non-zero 
background (see iCao et al.L [19971 §5.2 for a more detailed dis- 
cussion). 

Despite the improved resolution of the IGA over the orig- 
inal IRAS data, the resolution is not as good as that of MSX 
(18.3") and source contamination may still be an issue (see fig- 
ure|2|. Nearby sources in the IRAS point source catalogue (PSC) 
were identified in order to remove contaminating sources from 
the sky annulus, and flag sources where the flux measurement 
is potentially contaminated by other nearby sources. However, 
the PSC is not necessarily complete within the galactic plane 
and may not include cases where several sources cluster to give 
the appearance of a resolved source. Therefore, all sources were 
subsequently visually inspected and cases where the RMS target 
was dominated by other nearby sources were flagged as unus- 
able. 



2.3. MIPSGAL 

For sources in the inner galactic plane (5° < 1 <63° and 
298° < 1 <355°) to I b I < r, higher resolution data are available 
from the Spitzer Sp ace Telescope (SST ) MIPS Galactic plane 
(MIPSGAL) survey jCarev et all |2009|) . which has a resolution 
of 18" at 70 fim. Though fully reduced and calibrated mosaic im- 
ages have not yet been released for this waveband, observational 
data with basic reduction is available publicly for the whole sur- 
vey region. This was obtained and mosaiced in galactic coordi- 
nates with a pixel scale of 4 .8" pixel ' using the program mon- 
tage jBerriman et al. I l2006l) . Note that we do not consider the 
MIPSGAL 24 fim data as the majority of RMS sources are sat- 
urated, we already have MSX 21 /urn data and this waveband is 
not crucial to the determination of the total IR flux and therefore 
luminosity for our sources. 



The MIPS 70 jum detector suffers from a non-lin ear pixel re- 
sponse for surface brightnesses above ~66 MJy sr ' (iDale et akl 
2007) in wide-field mode or 1.5 mJy pixel \ corresponding to 
a point source flux of ~1 Jy. Though this is not a particularly 
large effect in extragalactic studies, the brightness of the galactic 
plane at these wavelengths results in background emission lev- 
els that are often of order 1000 MJy sr We therefore derive 
a pixel correcti on using th e following method, which is based 
in that used by iDale et akl (l2007i) . with the additional inclusion 
of errors on the data in our linear regression. For each standard 
star mea sured using aperture photometry in wide-field observa- 
tions by Gordon et al.l (120071, , their table 4), the measured flux 
in mJy can be calculated fro m the measured caUbratio n factor, 
the predicted flux (derived by lEngelbracht et al.L |2007|) and the 
mean calibration factor obtained from PSF fitting photometry 
of the same observations. The measured and predicted fluxes 
are then both converted to surface brightnesses by dividing by 
1.5 mJy pixeL', after which a linear bisector least-squares fit to 
those sources with fluxes > 1 Jy (see figure [3j is used to obtain 
the pixel correction of the form: 



(1) 



We find a = 0.558 + 0.051 and ^ = 1.11 + 0.02, which com- 



pares well with the results of iDale et al.l ( 2007 ). The differences 
between our correction and that o f lDaleetal.l(l2"007h come from 
the use of a linear bisector rather than simply a least squares 
linear fit. This correction was applied to all target images before 
photometry was undertaken. Though the MIPSGAL data used in 
this work has not un dergone the post-processing undertaken by 
iGordon et al.l (12007 ^. this primarily served to counter the slow 
detector response to the minimal background. RMS sources are 
generally much brighter and associated with strong background, 
so this issue is less of a concern. The systematic error introduced 
by the correction is usually of the order 15%, with the correction 
in the flux leading to an increase of ~ 62 % at 1000 Jy (see 
figure O. Though we are extrapolating the correc tion to fluxes 
well above those observed by iGordon et aTl (l2007l) . comparison 
of corrected MIPS 70//m and IRAS PSC 60/im fluxes (discussed 
in ^5. It returns a slope of 1.00 in log-space, suggesting that the 
correction is reasonable for all our sources. This is of a sim- 
lar order to that expected for such sources from more detailed 
corrections currently being undertaken by the MIPSGAL team 
(Paladini, private communication), and the error in this correc- 
tion dominates any other errors. 

As can be seen in figure [1] MIPSGAL 70 fim images show 
striping caused by differences in the background levels between 
observing scans across the sky. However RMS sources are usu- 
ally much brighter than the background and striping so the effect 
is not particularly strong on flux measurements. Therefore no 
attempts were made to remove it from the images. 



3. Two dimensional bacl^ground fitting photometry 

In order to perform background subtraction when the back- 
ground itself is spatially variable (e.g. see figure [U, a 2- 
dimensional fit to the background is the next logical progres- 
sion from a simple average. As in standard aperture photometry, 
an annulus around the source of interest is selected, though in 
this case the surface in the annulus is fitted with a 2-D second 
order polynomial using a least-squ ares multiple linear regres- 
sion method (see iBevingtoni 1 19691) . While a higher dimension 
fit could be used, it is questionable whether an annulus around a 
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Fig. 2. Example of a source for which nearby sources do not strongly contaminate the IGA flux (top, G015.0755-00.1212) and of a 
source which is too badly contaminated in the IGA for useful fluxes to be determined (bottom, G043. 1497+00.0272). Left: The IGA 
60 i-im images with 20 HIRES iterations. The standard aperture and sky annulus radii used in two-dimensional background fitting 
photometry are indicated using ellipses. The box indicates the area covered by the right hand image. Right: MSX E-band (21 fim) 
image with the IGA 60 fim aperture and sky annulus radii shown as before. For G043. 1497+00.0272 the RMS source is just visible 
on the edge of a large complex of much brighter sources which dominate the flux in the IGA aperture. 
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Fig. 3. Left: The measured and expected pixel surface brightness values derived from the data of iGordon et al.l (l2007h in order to 
obtain the MIPS 70 fim pixel surface brightness correction. The solid lines show linear bisector least-squares fit to fluxes > 1 Jy 
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source would uniquely constrain such a fit, or whether the varia- 
tion in diffuse background truly changes quickly on small spatial 
scales. Ideally, the background should vary on large scales com- 
pared with the size of a point source so that it is fully resolved. 

It is possible, however, that other nearby sources lie partially 
or wholly within the sky annulus (e.g. see figure|2l). In such cases 
the coordinates of these sources can be used to remove pixels up 
to a given radius around such sources from the area used for the 
fit, so that they do not affect it. It is not desirable to have too low 
an area from which to derive the fit, so a cut-off percentage is 
used below which no more sources are considered for removal. 

Having derived the fit to the pixels within the annulus, the fit 
is then removed from the original image and the total counts 
within the aperture calculated. The total error in the aperture 
counts is the combination in quadrature of the noise in the sky 
within the aperture, the error in the fit to the sky within the annu- 
lus and the errors due to the detector (e.g. read noise, gain etc.). 
With modern detectors it is often difficult to recover this last er- 
ror explicitly from the released mosaiced images, but it can be 
considered a limiting minimum error cr/„„, and can therefor be 
obtained through measurement of regions free from sources of 
emission. The total error is therefore given by: 



0-7-0/ - 



N„. 



Sfl„(Co - CfY 



Sfl„(Co - CfY 



+ a-t 



(2) 



where Co and c / are the pixel counts of the included pixels within 
the sky annulus of the original image and 2-D fit respectively 
while Nap and Nan are the number of pixels included in the aper- 
ture and annulus. 

In order to test this method, the results of aperture photome- 
try and 2-D background fitting photometry were cornp ared using 
TIMMI2 10.4 fim images from ' Mottram et a 042007) of sources 
in their 'clean' subsample (see figure H]). The images for these 
sources have approximately flat uniform backgrounds and iso- 
lated unresolved targets, so the results should be the same using 
either method. The same aperture and inner annulus radii were 
used as for the standard photometry for TIMMI2 observations 
(i.e. 12 and 15 pixels respectively) for both methods. However 
it was found that the 2-D aperture fitting photometry requires a 
larger sky annulus in order to constrain the fit properly. An outer 
radius of 30 pixels was therefore used for the aperture fitting 
photometry, as opposed to the outer radius of 20 pixels used for 
the standard aperture photometry measurements. The flux was 
also measured using the 2-D aperture fitting photometry at four 
locations on each image which were free of source emission and 
the mean calculated in order to obtain a limiting background flux 
for each source. The mean of this limiting background flux over 
all sources was then calculated to obtain the limiting error cr/,„,. 

The linear bisector least-squares fit to the photometric results 
for the subsample of TIMMI2 sources (see figure |4|i has a slope 
of 1.0 within the errors, and the correlation of aperture photom- 
etry to 2-D background fitting photometry fluxes is also 1 .0. The 
two methods are therefore consistent with each other within the 
calculated errors. In addition to comparisons between the two 
methods for the same image, comparisons were also performed 
between the original TIMMI2 fluxes and 2-D background fitting 
photometry obtained from images which had a diagonal ramp 
in the background counts and a large bright 2-D gaussian source 
near the target added to them. These returned very similar results 
to the first test, so the new method works well even in challeng- 
ing circumstances. 

In order to obtain accurate photometry for IGA data, which 
has an elliptical beam profile (see figure |2]l, both standard aper- 
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where m —1.01 with a —0.01 : 


and c =-0.015 with a =0.007 : 
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Fig. 4. Comparison of fluxes measured using aperture photome- 
try and 2-D background fitting photometry for the 'clean' sub - 
sample of TIMMI2 sources established in Mottr am et al. I (l2007h . 
The solid line shows a linear bisector least-squares fit to the data 
including errors while the dashed line indicates the line of equal- 
ity between the fluxes. The correlation between the data is given 
in the upper left corner of each plot. 

ture photometry and aperture fitting photometry routines were 
also developed to use elliptical instead of circular apertures and 
annuli. 



4. Results 

4.1. IGA Results 

Aperture fitting photometry was undertaken on all 1336 RMS 
candidate MYSOs identified as young sources, of which 602 
were too badly contaminated in IGA data for reliable results to 
be obtained, leaving 734 sources with acceptable measurements. 
The aperture and inner and outer sky annulus major axis radii of 
1 1, 28 and 46 pixels were used for the 60 fj.m data (as shown in 
figure 121) and 13, 30 and 50 pixels for the 100 jum data for all 
sources, after careful examination of isolated sources. The ratios 
of the major to minor axis radii used were 1.9 for 60 yum and 
1.3 for 100 jum. The limiting flux error cr/,„ was measured using 
background regions in 13 images for both 60 yum and 100 jum 
data, resulting in errors of 0.51 Jy and 5.97 Jy respectively. The 
much larger limiting error for the 100 fj.m images is probably due 
to the brightness of diffuse galactic dust emission at this wave- 
length. 

During their analysis, 'C ao et all ([1997) found that the IGA 
after 20 processing iterations overestimates the fluxes of iso- 
lated sources with respect to the IRAS PSC. The final flux re- 
sults therefo re include divisio n by the IGA(20)/IGA(1) factors 
measured bv lCaoet all (Il997h of 1.02 and 1.10 for IGA 60 jum 
and 100 fj.m data respectively. This correction results in a sys- 
tematic error in the fluxes due to the error in measurements of 
ICao et all (Il997h of 7% and 14% in the measured 60 fim and 
100 fj.m fluxes respectively. This systematic error in the fluxes is 
in addition to the error derived during the 2-D background fitting 
photometry, which is a random error Therefore these two errors 
are kept separate and the systematic error is not included in the 
results p resented in tab le [T] The correction factors were mea- 
sured bv lCao et al.l (1199 7) using relatively isolated sources with 
well defined backgrounds, so may well be larger for RMS can- 
didates which are often in complex star formation regions with 
high levels of background emission. 
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An example of these results is presented in table [T] while 
the full version of this table is available at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ The MSX PSC 
name for each target is given in column 1 and the IGA measure- 
ments and errors are presented in columns 2 and 3 respectively. 
The coordinates of the nearest IRAS PSC entry within 0.7' (the 
selection of which will be discussed in ^5.1b of the MSX coordi- 
nates of the target are given for comparison in columns 5 and 6. 
The offset between the MSX and IRAS PSC coordinates is given 
in column 7 while the IRAS PSC fluxes with errors are given 
in columns 8 and 9. The errors in the IGA fluxes were calcu- 
lated using equa tion^ and are sometimes quite small. However, 
ICao et all (1 19971) estimate that the overall photometric accuracy 
of the IRAS Galaxy Atlas is ~25%, so the dominant error in 
these results is probably from the image data rather than the mea- 
surements. 

Of the 734 targets with IGA fluxes, 682 have an IRAS PSC 
source within 0.7' but 172 of these PSC sources are upper limit 
detections at either 60 jim, 100 i^m or both. Therefore, IGA 
fluxes have been obtained for 224 sources which have no far- 
IR detection in the IRAS PSC at these wavelengths. Figure |5] 
shows the comparison between IGA and IRAS PSC fluxes for 
the 510 sources which have IGA fluxes and true IRAS PSC 
detections at 60 ;um and 100 fim. The IGA fluxes are, on av- 
erage, slightly larger than the IRAS PSC fluxes by a factor of 
-1.12 + 0.05 at 60 //m and -1.35 + 0.09 at 100 fim, however the 
fluxes in the IRAS PSC did not include a correction for hystere- 
sis at 60 yum, which when included tends t o increase t hese fl uxes 
and would account for the increase seen (ICaoetaLLfi^ . give 
IGA(1)/PSC = 1.12). The increase in the 100 fim data may be 
due to flux increases during HIRES processing as our sources are 
often in regions of considerable background emission which can 
cause interactions between this and the point source. However, 
the difference may also be due to more accurate background sub- 
traction. 

4.2. MIPSGAL Results 

2-D background fitting photometry was attempted for 843 Red 
MSX Source survey candidate MYSOs identified as young 
sources which lie within the MIPSGAL survey region. However, 
85 were either heavily saturated or have a much brighter source 
nearby so that no usable information could be obtained from the 
image (e.g. see the left hand plot of figure |6]l and 34 are visi- 
ble but saturated. Of the 724 sources where flux measurements 
were possible, 693 are detected and measurable in MIPSGAL 
images, though nearby sources had to be masked from the sky 
annulus for some of these sources (e.g. see the right hand plot of 
figure |6ll. In addition, 12 sources are not detected at 70 /im and 
therefore probably not MYSOs or H ii regions and 19 sources are 
non-detections due to there being bright sources or complexes 
nearby, and so have large upper limits on the flux associated with 
them. For an isolated unresolved source, aperture and sky radii 
of 13, 16 and 32 pixels were used, while these radii were set by 
inspection of the images for more complex sources. 

The point source function (PSF) of Spitzer MIPS observa- 
tions is circular but slightly non-gaussian and includes a low- 
level extended contribution. An aperture correction was there- 
fore calculated for the MIPSGAL images by performing pho- 
tometry on the image of an unresolved source with a relatively 
flat background using various aperture radii and inner and outer 
sky annulus radii, then dividing the results of photometry with 
a particularly large aperture which measures as close to the to- 



rt.UCI LU-iC 

Radius (") 


Sky Radius 
Inner (") Outer (") 


rt.UCI LU-iC 

Correction 


19.2 (4 pix) 
33.6 (7 pix) 
38.5 (8 pix) 
62.4 (13 pix) 


28.8 (6 pix) 
38.5 (8 pix) 
48.0 (10 pix) 
76.8 (16 pix) 


57.6 (12 pix) 
67.2 (14 pix) 
96.0 (20 pix) 
153.6 (32 pix) 


2.50 ± 0.42 
1.48 ± 0.17 
1.22 + 0.12 
1.14 + 0.11 



Table 2. Sample MIPSGAL 70 fim aperture correction values 
calculated using aperture fitting photometry for an isolated un- 
resolved source. 



tal source flux as possible. This ratio then provides an aperture 
calibration which is multiplied by the measured flux in order 
to arrive at the total flux for a given source. Calculation of the 
aperture correction was undertaken after pixel non-linearity cor- 
rection of the images (see ^2.3\ . A sample of correction val- 
ues are shown in table [2l and compare well to those obtained 
by ' Gordon et al.l (l2007h using model PSFs for blackbodies of 
10,000 K, 60 K and 10 K, particularly those for the 10 K black- 
body. The errors in the calculated corrections are generally be- 
low 10%, and are primarily due to the limiting error associated 
with the measurements used to derive the correction factors. 
The limiting error cr/„„ was found to be 2.68 Jy, measured as 
discussed previously, using 17 images with relatively flat back- 
grounds. 

For non-detections, four measurements of the background 
flux were taken using the standard aperture and annulus radii and 
the mean of the absolute values of the measurements calculated. 
A 3 cr upper limit on the flux for the source is then given by three 
times this mean. For those sources which are non-detections due 
to bright nearby emission, the flux was calculated in the same 
way as for a standard source, and this flux used as an upper limit 
to the true MIPSGAL 70 fim flux for the source. Though the 
saturation limits of MIPSGAL 70 fim images have yet to be fully 
characterised, it is certainly a problem for bright RMS sources 
or sources in or near bright regions of the galactic plane. 

Example MIPSGAL 70 fim 2-D background fitting pho- 
tometry results are shown in column 4 of table [1] while 
the full version of this table is available at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or 
via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/, All sources 
which either have a flux measurement or are non-detections for 
whatever reason are included in the full table. 

While the resolution of SST MIPS 70 fim observations is 
much better than previous surveys (18"), allowing direct com- 
parison with the MSX su rvey, higher resolution mid-IR imaging 
of RMS sources (e.g. see lMottram et aUl2007l) have shown that 
source confusion may still be an issue. Every effort has been 
made to ensure that only the source of interest is included within 
these flux measurements, but it is possible that some unresolved 
sources may in fact contain more than one emitting object. 

5. Discussion 

5.1. Comparison of IRAS and MIPSGAL photometry 

In order to better evaluate the MIPSGAL results, a search for 
IRAS PSC entries associated with RMS sources was conducted. 
For each source with a MIPSGAL 70 fim detection, the near- 
est IRAS PSC entry with detections at both 60yum and 100 fim to 
the MSX coordinates of the RMS source was identified. Figure|7] 
shows the mean IRAS 60 ftm to MIPSGAL 70 fim flux ratio as a 
function of offset between the IRAS and MSX coordinates. From 
this figure, as well as visual inspection of MSX and MIPSGAL 
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Table 1. Example of measured IGA 60 //m, IGA 100 yum and MIPSGAL 70 /im photometry, see text for details. A 
full version of this table is available online at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via 
|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/, 



Object IGA 60yum IGA lOOyum MIPSGAL 70yum Nearest IRAS PSC Entry IRAS 60//m IRAS lOOyum 
Flux (ly) Flux (ly) Flux (Jy) R A (12000) Dec (12000) Offset (') Flux (ly) Flux (Jy) 



GO 11.0732-00.2295 
GOl 1.1 109-00.4001 
GOl 1.1723-00.0656 
GOl 1.9454-00.0373 
G012.5932-00.5708 
GO 12.9090-00.2607 
G013.2097-00.1436 
G013.6562-00.5997 
G013.8885-00.4760 
G014.0329-00.5155 
G014.2071-00.1105 
G014.2166-00.6344 
G014.6703-00.4795 
G015.0755-00.1212 
G015. 0939+00. 1913 
G015.4805-00.3900 
G016. 1438+00.0074 



134.48±2.14 
1145.32±1.86 
181.62±2.85 
1423.15±3.95 
94.24±1.75 
3212.46±15.41 
1214.14±4.21 
403.97±2.27 
661.08±4.96 
265.45+16.31 
346.53+4.85 
70.15+4.50 
138.96+1.63 
249.51+4.06 
127.00+2.36 
175.60+1.20 
44.23+4.26 



132.32+12.52 
2323.30+10.73 
299.14+10.16 
1884.77+14.14 
229.58+9.25 
8116.61+27.53 
2999.06+14.72 
989.68+9.72 
861.81+26.79 
1195.19+8.10 
566.17+14.71 
398.38+13.86 
594.71+35.98 
575.18+14.22 
480.57+8.71 
712.90+6.77 
161.02+14.41 



59.35+3.20 
888.59+3.56 
114.20+3.16 
1119.95+3.27 

17.77+4.01 
2136.63+5.53 
1118.61+3.59 
403.22+3.12 
169.60+3.57 

70.79+3.64 
201.11+3.48 

57.22+3.21 

40.19+3.47 
211.48+3.29 

75.13+3.14 
<38.80 

85.67+3.17 



: 10:48.7 
:11:33.2 
:10:25.2 
:11:53.0 
:15:10.0 
: 14:39.0 
: 14:47.8 
: 17:24.4 
: 17:24.7 
:17:50.7 
t 

18:18:38.5 

t 

18:18:27.6 

t 
t 

18:20:03.7 



-19:27:42 
-19:30:39 
-19:17:51 
-18:36:20 
-18:17:30 
-17:52:03 
-17:32:48 
-17:22:13 
-17:06:26 
-16:59:55 

-16:53:33 

-15:53:36 



0.210 
0.219 
0.078 
0.116 
0.105 
0.183 
0.511 
0.027 
0.023 
0.040 

0.088 

0.086 



-14:53:31 0.260 



101.00+12% 
853.00+16% 
73.30+13% 
1180.00+17% 

<57.90 
2230.00+23% 
877.00+22% 
422.00+21% 
262.00+16% 
< 145.00 

34.40+18% 

181.00+20% 

98.20+9% 



294.004 
2090.00: 

<330, 

1710.00; 

<624, 
6310.00: 
1850.00: 

1 140.00: 

407.00i 
923.004 



:18% 

tl7% 
00 

tl5% 
00 

t22% 
t27% 
tl7% 

:36% 
:25% 



<2820.00 
<477.00 

233.00+17% 



No IRAS PSC entry within 0.7' of MSX target coordinates. 
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Fit is of the form 








LOG(Y) = m.LOG(X) + c 






where m =0.98 with o =0.01 






ond 


c --0.05 with a -0.02 





100 10^ 
IGA 60/im Flux (Jy) 



10" 




100 10^ 
IGA 100,am Flux (Jy) 



Fig. 5. A comparison of measured IGA 60 jum fluxes with IRAS 60 //m (left) and IGA 100 /zm fluxes with IRAS 100 jum (right) PSC 
fluxes for 510 sources with PSC detections. The solid line shows the linear-bisector least-squares fit to the data while the dashed 
line indicates the line of concordance. 



images, the radius out to which the nearest IRAS PSC candidate 
is considered a probable counterpart to the RMS candidate was 
set at 0.7'. This is more than twice the resolution of MSX (18.3") 
and almost three times the quoted IRAS PSC cross-scan posi- 
tional error (15.6")- This limit was also used to identify probable 
IRAS counterparts for sources with IGA observations (see ^4.11 ). 

Figure [8] shows the comparison between the measured 
MIPSGAL 70 jum photometry and IRAS PSC 60 yum (left) 
and 100 yum (right) fluxes for all 242 sources with measurable 
MIPSGAL fluxes and detections at both 60 //m and 100 yim in 
the IRAS PSC within 0.7'. Sources which have upper limit IRAS 
PSC fluxes (i.e. cr - 100%) at either 60 yum or 100 ycca are not 
included. Linear bisector least-squares fits to the data give rela- 
tionships between the MIPSGAL 70 yum and IRAS fluxes of the 
form shown in equation [1] with powers equal to one within the 
errors and flux ratios of F/r^s 60 / Fm/ps vo ~ 1.32+0.21 
and F/j?As loo / '^mips 70 ~ 2.82+0.65. The expected 
flux ratios assuming a blackbody of temperature 50K 
(or 40K) are F/ras 60 / F^/ra 70 = 1-12 (or 0.95) and 
F/RAS 100 / Fm/ps 70 - 1-88 (or 2.43), so the observed ratios 



may be in part due to the differences in filter profiles between 
IRAS and MIPS. However, the mean ratios derived from fits 
to the SEP models o f RMS sources using the model fitter of 
iRobitaille et all (l2007h are 0.93 and 1.07 for F/ras 60 / ^mips 70 
and F/RAs loo / Fm/ps 70 respectively (Mottram et al., in prep.). 
The model fitter convolves the models with the observed filters 
using the standard colour correction, so the model fluxes and 
flux ratios are directly comparable to the observations. Though 
the IRAS filters are much wider than the MIPS filter, this effect 
is included in the ratios produced by the model fitter as well, 
so the large ratio between IRAS PSC and MIPSGAL fluxes 
is probably due to contamination in the larger IRAS beam, 
particularly in the 100 yum band. The smaller model ratio of 
F/pAS 60 / F/RAs 100 than would be expected for a blackbody of 
relevant temperature is probably the result of the distribution of 
dust temperatures within the SED models. 

The general spread of sources in figure [8] is due to several 
factors. The SED of all young sources will be similar, but not 
identical, leading to some fluctuation of the flux ratios. As men- 
tioned above, sources in more crowded regions are likely to have 
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6 4 2 -2 -4 -6 6 4 2 -2 -4 -6 

Arcminutes Arcminutes 

Fig. 6. Examples of MIPSGAL 70 fim images for RMS targets. The apertures and sky annuli used in photometry are shown with 
open circles while excluded regions are indicated by line-filled circles. Unphysical pixels are shaded orange in the online version of 
this paper. Right: a source (G059. 7856+00. 0734) where nearby bright sources would confuse any flux measurement. Left: a source 
(GOIO. 6207-00. 3 199) with brighter nearby sources which would be confused in IGA observations (60 //m aperture shown by an 
ellipse). 



40 r 



30 : 



100 



20 - 



10 
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0.5 1.0 
Offset (') 
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Fig. 7. The mean ratio of IRAS PSC 60 jum flux to measured 
MIPSGAL 70 yum flux for the nearest IRAS PSC entry to the 
target MSX coordinates as a function of off'set. The solid line 
shows the linear bisector least-squares fit to the data below 0.7', 
the dotted line shows the fit to the data between 0.7' and 1.5' and 
the dashed line shows the fit between 0' and 1.5'. The vertical 
dot-dashed line indicates 0.7', above which offset the correlation 
between RMS targets and the IRAS PSC becomes weaker. 

more confused IRAS fluxes, leading to both larger flux ratios 
and larger variation. For sources closer to the galactic centre, the 
line of sight of observations passes through more of the galactic 
plane, so there is more chance of confusion. 



+ 




1 10 100 

'^MSX 21 / '^MSX 8 

Fig. 9. Mid and far-IR colour-colour plot using MSX and 
MIPSGAL 70 yum fluxes. Sources with RMS identifications of 
'YSO' are shown with '+' signs, 'Hii / YSO' with squares, 'Hii 
region' with filled circles and 'Young / Old star' with triangles. 
The vertical solid line indicates the RMS Fmsx i\ I Fmsx 8 se- 
lection cut, while downward arrows indicate sources with non- 
detections at 70 yum. Note that sources with fluxes above this cut 
within the 1 cr errors are still included in the sample. 



5.2. Comparison of MSX and MIPSGAL colours 

Figure |9] shows a colour-colour plot using MIPSGAL and MSX 
fluxes for all 685 sources which have MIPS 70 yum detec- 
tions and Fmsx 8 > in terms of RMS identification (see 
92.11 ). While young sources have both cold and warm dust as- 
sociated with them, more evolved sources such as AGB stars 
will only be associated with strong emission from warm dust 



(see e.g. Ivan der Veen & Habind. [19881) . The five sources with 
Fmips 70 / Fmsx 21 < 1 are therefore unlikely to be MYSOs. 
Of these, two (G029.4332-h00.1542 and G03 1.9844-00.4849) 
are bright but relatively isolated in GLIMPSE images and so 
are most likely evolved stars. For comparison the SEDs of 
G029.4332+00.1542 and G034.7569+00.0247, a typical YSO, 
are shown in figure [TO] It can be seen that inclusion of far- 
IR flux data is important for discriminating between young 
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Fit is of the form 
LOG(Y) = m.LOG(X) + c 
where m =1.00 with a =0.03 
ond c =0.1 2 with a =0.07 



100 10^ 
MIPSGAL 70/Lim Flux (Jy) 



10" 




Fit is of the form 
LOG(Y) = m.LOG(X) + c 
where m =0.99 with a =0.06 
ond c =0.45 with a =0.10 



10 



100 10' 
MIPSGAL 70/im Flux (Jy) 



Fig. 8. A comparison of measured MIPSGAL 70 fim fluxes with IRAS 60 yum (left) and 100 fim (right) PSC fluxes for 243 sources 
with IRAS PSC detections within 0.7' of the MSX target. The solid Une shows the linear-bisector least-squares fit to the data while 
the dashed line indicates the line of concordance. 




1000 



G034.7569+00.0247 



10 



100 



A (/.im) 



1000 



Fig. 10. The SEDs of G029.4332-i-00.1542, an evolved star with similar mid-IR colours to young sources and the 'well-behaved' 
YSO G034.7569-h00.0247. The filled circles indicate detected fluxes while the trian gles indicate upper limits. The data used 
include 2MASS MSX IRAS, SCUBA s ub-millimetre and MAMBO 1.2mm fluxes (I Skrutskie et"aU 120061: lEgan et all 120031: 
iBeichman et al.L Il988t iDi Francesco et aL . 2008; Beuth eret all l2002bh in addition to those presented in this paper. These and 
the SEDs of all young RMS sources will be discussed further in Mottram et al. in prep. 



and evolved sources. A further two (G034.8625-00.0629 and 
G334. 7202 -1-00.1762) appear similar to the previous sources in 
GLIMPSE images but are also associated with unresolved ra- 
dio detections, so are most likely PNe. All four sources show 
^^CO detections, probably due to diffuse clouds either behind 
or in front of them. The last (G3 16.7754-00.0447) is certainly 
in a star forming region as both MSX and GLIMPSE images 
show multiple sources and diffuse emission nearby. However 
no point sources are evident in any of the Spitzer IRAC bands, 
so this source is either a more evolved Hii region or a photo- 
dissociation region (PDR) which has bright mid-IR emission. 
RMS sources already identified as evolved are not included in 
this paper, and have not had their MIPSGAL 70 /urn fluxes mea- 
sured. 

In general the H ii regions have larger Fmsx i\ I Fmsx 8 and 
Em/PS 70 / Fmsx 21 ratios than the YSOs, which is consistent with 
the H II regions being more embedded than the YSOs. However, 
there is a population of YSOs which have similar colours to the 
Hii regions. 



6. Summary and conclusions 

In this paper, we have presented a new technique for background 
subtraction in aperture photometry using a 2-D fit to the back- 
ground within the sky annulus. This was then used to obtain pho- 
tometry from MIPSGAL 70 yum mosaic images for 724 young 
RMS sources and for 734 such sources using IRAS Galaxy Atlas 
60 iim and 100 yum images. Overall, far-infrared fluxes have been 
obtained for 1113 of the 1336 candidates examined, of which 
370 have no IRAS PSC entry and a further 374 have only upper 
Umits in the IRAS PSC at 60 //m and/or 100 //m. In addition, 4 
objects have been identified as evolved due to their MIPS fluxes. 
This is one of the first use s of MIPSGAL 70 [i m data (another 
example being the work of IChapin et al.l (l2008h ) and one of the 
largest sets of measurements of far-IR sources in the galactic 
plane. Comparisons with the IRAS PSC, where applicable, show 
reasonable agreement, though the MIPSGAL fluxes are strictly 
superior due to improved resolution, and the IGA data includes 
a hysteresis correction at 60 yum which the IRAS PSC did not. 
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These far-IR flux measurements, along wifli other photomet- 
ric data, provide information about the spectral energy distri- 
bution (SED) of young RMS sources, and thus allow the de- 
termination of the luminosities of this sample of sources. This 
will allow the systematic study of the properties of a large, well- 
selected sample of MYSOs as a function of luminosity, and will 
be the subject of forthcoming papers. Information on all sources 
discussed in this paper is available via the RMS database, which 
can be found at www.ast.leeds.ac.uk/RMS/. 
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